Matched-filtering and parameter estimation of ringdown waveforms 
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' Using recent results from numerical relativity simulations of non-spinning binary black hole merg- 

ers we revisit the problem of detecting ringdown waveforms and of estimating the source parameters, 
considering both LISA and Earth-based interferometers. We find that Advanced LIGO and EGO 
' could detect intermediate-mass black holes of mass up to ~ 10'^ Mq out to a luminosity distance 

I of a few Gpc. For typical multipolar energy distributions, we show that the single-mode ringdown 

QQ , templates presently used for ringdown searches in the LIGO data stream can produce a significant 

event loss (> 10% for all detectors in a large interval of black hole masses) and very large parameter 
estimation errors on the black hole's mass and spin. We estimate that more than ~ 10^ templates 
would be needed for a single-stage multi-mode search. Therefore, we recommend a "two stage" 
search to save on computational costs: single-mode templates can be used for detection, but multi- 
mode templates or Prony methods should be used to estimate parameters once a detection has 
5h , been made. We update estimates of the critical signal-to-noise ratio required to test the hypothesis 

that two or more modes are present in the signal and to resolve their frequencies, showing that 
second-generation Earth-based detectors and LISA have the potential to perform no-hair tests. 
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I. INTRODUCTION 

O .... . n, 

' Astronomical observations provide us with a large number of black hole candidates [![. Stellar mass candidates 
, (with mass M ~ 5 — 20Mq) are found in X-ray binaries in our galaxy. Since the discovery of quasars and other active 
galactic nuclei, we have growing evidence that supermassive black holes (SMBHs, with M ^ 10^ — 10^ Af©) should 
reside in the cores of almost all galaxies, including our own. There is also mounting observational support to the idea 
that intermediate-mass black holes (IMBHs) may fill the gap between these two classes of astrophysical objects. 
?H I It is commonly believed that the ultimate test of the "black hole hypothesis" will come from gravitational wave 
observations (see eg. Lasota's entertaining account of exotic alternatives to astrophysical black holes 0, and Hughes' 
lectures Q for a relativist's perspective on present and future observational tests of the black hole hypothesis). 
From the point of view of an astrophysicist, black holes are not particularly interesting: "black hole candidates" are 
characterized only by their mass and angular momentum, and no compact object with mass M > 3Mq has shown 
any feature that would allow us to attribute to it any other property other than mass and rotation. For a relativist, 
black holes (being vacuum solutions of the field equations) are much more exciting: they are unique, "clean" probes 
of the structure of spacetime in strong-gravity conditions. 

An important identifying dynamical feature of a black hole are its characteristic damped oscillation modes, called 
quasinormal modes (QNMs) j^]. Any compact binary merger leaving behind a black hole produces a gravitational 
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wave signal that, after the black hole's formation, can be decomposed as a superposition of exponentially damped 
sinusoids. By analogy with the ordinary vibrations of a bell, this signal is known as "quasinormal ringing" or ringdown. 

The fact that all information is radiated away in the process leading to black hole formation, so that astrophysical 
black holes in Einstein's theory are characterized completely by their mass and angular momentum, is known as the 
"no-hair theorem" . For this reason the ringdown signal is very simple: mass and angular momentum are enough to 
characterize the black hole's oscillation spectrum (see Ref. [5] for details), and the detection of ringdown waves may 
allow us to identify a black hole and determine its properties. A measurement of the frequency and damping time of a 
single QNM can be used, at least in principle (see below), to determine both the mass and angular momentum of the 
hole @,S0, 9- Detection of more than one QNM would allow us to test consistency with the black hole hypothesis 
- or in more colloquial speaking, to "test the no-hair theorem" [1,0. 

The detectability of ringdown waves, as well as their use to measure black hole properties and test relativity, 
depend mainly on two factors: (i) the fraction of the black hole's mass radiated in ringdown waves (i.e. the "ringdown 
efficiency" Crd), and (ii) the detector's sensitivity in the frequency band of interest. In Ref. some of us carried out 
a detailed study of ringdown detectability and no-hair tests with the planned space-based Laser Interferometer Space 
Antenna (LISA). The ringdown efhciency Crd, the multipolar energy distribution of the radiation and the dimensionless 
angular momentum j of the final black hole were considered as free parameters, to be varied within a certain reasonable 
range. We confirmed and refined earlier estimates by Flanagan and Hughes [8l], showing that the ringdown signal-to- 
noise ratio (SNR) is usually larger than the inspiral SNR for the typical SMBH masses {M > 10^ ^0) inferred from 
astronomical observations of nearby galaxies. 

Binary black hole simulations are now carried out by different groups all over the world. Large-scale simulations of 
non-spinning unequal-mass binary black hole mergers were recently used to provide reliable estimates of e^d, of the 
final angular momentum j and of the multipolar energy distribution [iq| . Ref. jT(l| also showed that two or more modes 
are excited to comparable amplitudes in a binary black hole merger whenever the binary's mass ratio q = 7712/7711 1. 
The significant excitation of different multipolar components can be used to extract useful information on the geometry 
of the ringing object, and to perform no-hair tests (in the sense explained above). 

In this paper we revisit the analysis of Ref. fs^ taking into account these recent results from numerical relativity, 
and we extend that study to include planned and presently operating Earth-based detectors (LIGO, Virgo, Advanced 
LIGO and EGO). Many authors have recently stressed t he p otential of Earth-based detectors for measuring ringdown 
waves in the IMBH mass range For example, Ref. fi3\ matched an equal-mass 

merger waveform from numerical relativity to a Post-Newtonian inspiral and showed that the resulting SNR for a 
single LIGO detector peaks at M ~ 150 Mq, which is well within the IMBH mass range. The results presented 
in this paper confirm that ringdown waves may be used to provide conclusive evidence for the existence of IMBHs 
and (even more importantly) to accurately measure their parameters. Initial LIGO (Virgo) may detect IMBHs with 
M < 400 Mq {M < 800 Mq, respectively) out to a luminosity distance Dl ~ 100 Mpc. Advanced Earth-based 
interferometers would extend the luminosity distance at which ringdown events are detectable by a factor 10 for 
Advanced LIGO, and ~ 100 for EGO (see Fig. [2] below). Advanced LIGO and EGO could detect IMBH ringdowns 
out to cosmological distances with large SNR, and allow precision measurements of the IMBHs' properties. 

Higher multipolar components of the radiation can be significantly excited (IQ] , and this calls for a critical revision 
of the matched-filtering techniques used for ringdown searches. Present ringdown searches in the LIGO data stream 
use templates consisting of only one QNM (see also Ref. f23!| for data analysis methods to resolve neutron star 
ringdowns from instrumental glitches). 

In this paper we address this problem and give a preliminary answer to the following questions: 

(1) Given estimates of the multipolar energy distribution in non-spinning binary black hole mergers, how many 
events would we miss in a search with single-mode ringdown templates? The answer is quantified in Fig. [Jl 
where we compute Owen's "minimal match" [IJ] as a function of the black hole mass AIq measured in the 
source frame. The event loss is larger than 10% whenever the minimal match is larger than 0.035 (details are 
provided in the body of the paper) . From the figure we see that this happens in a significant mass range for 
all Earth-based detectors. We will show below that similar conclusions apply also to the planned space-based 
interferometer LISA, although in a completely different mass range. Furthermore, we will show that single-mode 
templates can produce a large bias in the estimation of the black hole's mass and spin. Our conclusions should 
be rather conservative, because from perturbation theory and from numerical relativity results we expect higher 
multipoles to be more excited for black hole binaries with spin and large mass ratios. 

(2) How many templates would we need for searches using two- mode templates? We estimate that, as compared 
with single-mode searches, the number of templates needed for a two-mode search would increase by roughly 
three orders of magnitude. A more detailed data analysis study (eg. using better template placing techniques, 
along the lines of [25l|) could be useful to reduce computational requirements. 



3 




FIG. 1; Owen's minimal match, as defined below Eq. p.5|l . for different Eartli-based detectors. For tliis illustrative calculation 
we assume that the Kerr parameter of the final black hole is j = 0.6, and that the relative amplitude of the second mode is 
A — 0.3. We also set the phases in Eq. (|3.ip to be 0i = (/>2 = (thick lines) or = 0, 02 = tt (thin lines). The black circle 
and the red square mark two cases that we study in more detail below: a Mo = 100 Mq black hole as observed by LIGO and 
a Mo — 200 Mq black hole as observed by Advanced LIGO, respectively. 



(3) How strong must a ringdown event be if we want to perform no-hair tests or resolve nonlinear contributions [26[ 
to the ringdown waveform? More precisely: what is the minimum SNR needed to resolve QNMs? Below we 
address different aspects of this problem (frequency, damping time and amplitude resolvability require different 
SNRs). Our results suggest that prospects to resolve QNMs are quite optimistic for both LISA and second- 
generation Earth-based detectors. 

The plan of the paper is as follows. In Sec. |lT] we discuss the potential of Earth-based detectors to measure 
ringdown waves from solar-mass black holes and IMBHs. In Sec. IIIII we look at the event loss and bias in parameter 
estimation induced by searching for a two-mode waveform with a single-mode template. In Sec. IIVI we revise our 
previous estimates [5] of the critical SNR required to perform no-hair tests. Sec. fVl contains a brief summary of our 
conclusions. Some technical details are presented in the Appendices. In Appendix [X] we estimate the number of 
templates required to detect multi-mode waveforms, and in Appendix[B]we provide a (somewhat optimistic) estimate 
of the critical SNR required to test the hypothesis that two modes are present in a ringdown signal. 



II. RINGDOWN DETECTABILITY BY EARTH-BASED DETECTORS 



In this Section we study the detectability of ringdown waves by present and planned Earth-based interferometers: 
LIGO, Virgo, Advanced LIGO and EGO. Our analysis complements that of Ref. where a similar study was per- 
formed for the planned space-based interferometer LISA. Detectability can be assessed by computing the (maximum) 
SNR p achievable by matched-filtering: 



p^,/ih\h), (2.1) 
where the scalar product between two waveforms is defined as 0, [13] 





Here 

f -t-oo 



h{f)s / e"'-h{t)dt, (2.3) 

is the Fourier transform of the waveform h{t)^ and Sh{f ) is the noise power spectral density (PSD) of the detector. For 
the initial LIGO and Virgo PSD we use the analytic approximation of Ref. [2g|. We assume that the PSD Sh{f) = oo 
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for / < /s, where the low-frequency sensitivity cutoff fs = 40 Hz (for initial LIGO) and fs = 20 Hz (for Virgo). 
For Advanced LIGO we consider the broadband configuration PSD given in Ref. and for the projected PSD of 
EGO we follow Appendix C of Ref. [30] ■ For LISA, we model the PSD (including white-dwarf confusion noise) by the 
semianalytic approximation used in [a. |3]| . with a low- frequency cutoff /s — 3 x 10~^ Hz. 

An estimate of the maximum detectable black hole mass for LISA and Earth-based detectors can be obtained by 
equating the fundamental I ^ m ~ 2 QNM frequency (as tabulated and fitted in Ref. Q) to /s. Since QNM frequencies 
scale as 1/M, this sets a limit on the maximum detectable black hole mass. The results of this estimate, which are 
mildly dependent on the black hole's rotation parameter j, are listed in Table IH The mass range accessible to Earth- 
based interferometers can increase dramatically with relatively small (but technically difficult) improvements in the 
low-frequency sensitivity threshold. For example, a detector with /s=10 Hz can detect IMBHs with M > 10^ Mq 
even if they are non-spinning. For the typical spins that should result from a binary black hole merger (j ^ 0.6), 
a low-frequency sensitivity threshold at 40 Hz means that initial LIGO will not be sensitive to QNM ringing from 
IMBHs with mass > 400 Mq. The detectable mass range increases as the spin of the remnant black hole gets larger. 



TABLE I: Upper limit on the detectable black hole mass for which the fundamental QNM with I = m = 2 would be detectable 
by LISA and Earth-based interferometers, for different dimensionless spin parameters j. A reasonable estimate of fs is ~ 10 Hz 
for EGO, and ~ 20 Hz for Advanced LIGO (see Ref. :305). 
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FIG. 2: Left: ringdown SNR for LIGO, Advanced LIGO, Virgo and EGO at 100 Mpc. Thin lines refer to the EF convention, 
thick lines to the FH convention (see text). Right: luminosity distance (in Mpc) for detectability of the ringdown signal with 
a SNR of 10 (for clarity, here we only show results for the FH convention). 



In the left panel of Fig. [2] we show the typical sky-averaged SNR of Earth-based detectors for black hole ringdowns at 
luminosity distance Dl = 100 Mpc as a function of the black hole mass in the source frame Mq (related to the mass M 
in the detector frame by Af = (1 -I- z)M[), where z is the cosmological redshift). The actual value of the ringdown SNR 
changes slightly dependingon the way we compute the Fourier transform in (|2.3p . There are two common conventions 
in the literature (see Ref. for more details). The Echev err ia- Finn (EF) convention (in the terminology of Ref. 
is shown by thin lines in the left panel. It assumes that the ringdown waveform h[t) ~ e"*/"^ sin(27r/i) is zero before 
some starting time (say t ~ Q). Alternatively, we can adopt the Flanagan- Hughes (FH) doubling prescription (thick 
lines in the figure): we assume that the waveform for t < is identical to that for i > except for the sign oi t/r 
(i.e., we replace e"'/"^ by e"'*'/"^) and we divide the resulting SNR by to compensate for the doubling Q. The 
left panel of Fig. [2] shows that adopting one or the other convention does not change the results significantly. In the 
remainder of the paper, we will adopt the EF convention, since it is generally more straightforward to implement in 
matched-filtering applications. 
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The right panel of Fig. [2] shows the luminosity distance out to which ringdown would be detectable by each detector 
with a SNR of 10 (assuming a flat, A-dominated cosmology, in accordance with the latest observational data). In the 
Figure we use typical values from simulations of non-spinning binary black hole mergers [13] , assuming that the final 
black hole has angular momentum j — 0.6 and that the ringdown efficiency is ej-d ~ 3%. The results would not change 
much had we made different assumptions. In particular, the SNR scales with efficiency as p ^ ^erd- If the merging 
black hole binary has mass ratio q = •012/1711 < 1, the ringdown SNR would decrease (to a very good approximation) 
by a factor Aq/ (1 + q)^ compared to the equal-mass case [lOj. Fig. [5] can be compared with Figs. 7 and 8 in Ref. 
showing the inspiral and ringdown SNR of LISA at different values of the cosmological redshift (for clarity, in Fig. [5] 
we do not show the inspiral SNR). We also note that Figure 2 assumes matched filtering with one mode only, and 
that (following results from numerical relativity simulations) we assumed this mode to be the fundamental mode with 
I = m = 2. Inclusion of higher multipoles, the main topic of this paper, would only increase the sky-averaged SNR by 
small factors of order unity, but (as we will show) it may have important implications for a matched-filtering detection 
of the signal. 

The left panel of Fig. [2] is in good agreement with the results for LIGO and Advanced LIGO shown in Fig. 2 of 
Ref. Our results for Advanced LIGO are slightly different because we use a more accurate model for the PSD. 
For Virgo, our results agree with those in Ref. J^. The right panel clarifies the potential of present and future 
interferometric detectors to detect ringdown waves. At the present sensitivity, LIGO could detect waves from IMBHs 
of mass M ~ 10^ M© out to a distance of about 100 Mpc. Virgo spans roughly the same distance range, but it 
should be able to detect larger IMBH masses when the low- frequency design sensitivity is met. Advanced LIGO 
will extend the distance range out to a few Gpc (a factor ^ 10) and it should be able to detect IMBHs of mass 
10^ Mq < Mq < 10^ Mq. ego would be a remarkable tool to detect IMBHs of mass as large as ~ 10^ Mq out to 
distances of ^ 10 Gpc and larger, and to measure their properties with remarkable precision. 



III. EVENT LOSS AND BIAS IN PARAMETER ESTIMATION USING SINGLE-MODE TEMPLATES 

A common choice to search for signals of known form in noisy data is matched-filtering. Matched-filtering works by 
cross-correlating the signal with a set of theoretical templates^ . Current searches for ringdown signals in the LIGO 
data stream [22| and in resonant bar detectors ^33j] make use of the simplest theoretical model: a single damped 
sinusoid. Here we try to assess the performance of such a template to detect a superposition of damped sinusoids. 
More precisely, we estimate the number of events we would miss by using single-mode templates to detect multi-mode 
signals, and we try to quantify the bias in measured parameters induced by the use of such templates. 

The optimal way of addressing the performance of single-mode ringdown templates would be to use our present best 
guess for the "true" waveform emitted by a binary black hole merger: the wavetrain obtained by stitching together 
the Post-Newtonian predictions for the inspiral and the best available numerical relativity waveforms. Finding the 
optimal wayto perform this stitching is in itself a difficult problem, now actively pursued by many research groups 
[lol [TtI . [1^ [3^ . Our purpose here is to make some general points of principle, exploring semi-quantitatively the typical 
event loss and parameter bias induced by single-mode templates. For this reason, and to keep the analysis as simple 
and general as possible, our "true" waveform will be taken to consist of a two-mode QNM superposition. Therefore 
we assume that the following two-mode signal is impinging on the detector, starting at time t = 0: 

hit) = ylie"(''^i/'3i)'sin(27r/it-0i)-f yl2e~^''-''"/'5"^*sin(27r/2i-(/)2) (3.1a) 
= Ai Te-f'^^i/'^i)* sin (27r/ii - 0i) -t- yle"^'^-^^^/'^^)* sin (27r/2i - h)] ■ (3.1b) 

Here (z = 1, 2) denotes the oscillation frequency of each QNM, Qi = irfiTi is the quality factor of the oscillation, 
and Ai is the oscillation amplitude of mode i. For reasons that will become apparent in the following, in the second 
line we found it convenient to express the signal in terms of the relative amplitude of the two modes A = A2/A1. 
The Fourier transform of Eq. (j3.ip in the EF convention reads: 



^iJj=^i „rf2 AAffn — 7T( — f2 1 f2\^2i ^ (1 ^ ^) • (3.2) 



Qi [2/iQi COS cf)i - (/i - 2ifQi) sin ( 
7r[/f-4z//iQi+4(-/2-f /2)Q2] 

Matched-filtering works by cross-correlating the detector's output with a set of templates. Consider a one-mode 
bank of templates such as those presently used in ringdown searches [22| : 

T{\} = e-''f^^'^^^sin[2T:fTit-ta)~(j)T], iov t > to . (3.3) 



^ For alternative techniques especially designed to extract damped sinusoidal signals from noise, see Ref. [3^ and references therein. 
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where A is an (iV-dimensional) vector of template parameters, and T{\} = for < < to- In this particular case, 
A = (/t , Qt ,</'t) and = 3 (in principle we could use to as an additional parameter, but we verified by explicit 
calculations that setting to — doesn't sensibly affect the fitting factor and the estimated values of the other param- 
eters). The Fourier transform of this template is similar to those in Eq. (13. 2|) . with an additional exponential factor 
depending on the starting time to: 



~ ^ Qt [^/tQt cos./'t - (/r - 2^/QT) sin^r] 

7r(/|-4z//TQT + 4(-/2 + /|)Q|) '''P 



MfT - 2ifQT)to 



(3.4) 



The performance of the template (|3.3p to detect the two-mode waveform (|3.ip can be computed with the help of 
the fitting factor (FF) first introduced by Apostolatos [27|: 

FF = max^JM^. (3.5) 

A 



^iT{\}\T{X}){h\h) 



where the scalar product between two waveforms has been defined in Eq. (|2.2p . An equivalent quantity often used 
in the literature is Owen's "minimal match", defined as (1 — FF) [l^]. From the definition, Eq. (|3.5|) . it is clear that 
an overall normalization constant (say, Ai) does not affect calculations of the FF: this is the reason for introducing 
the relative amplitude A defined in (|3.ip . The FF measures the degradation of the SNR due to cross correlating an 
arbitrary signal h{t), such as the two-mode signal p.ip . with all filters T{X} in the template bank. The effective SNR 
that can be obtained by matched-filtering is 

PMF = max ^ = FF X p , (3.6) 



A 



(T{A}|T{A}) 



where p has been defined in Eq. (|2.ip . For gravitational wave detection the SNR is proportional to the inverse of the 
luminosity distance to the source, while the event rate (scaling with the accessible volume) is proportional to the cube 
of this distance. Therefore, given the FF, we can compute the event loss as 

event loss = f - FF^ (3.7) 

For detection purposes one usually requires FF > 0.965, corresponding to a loss of less than 10% of the events that 
could be potentially be detected by a "perfect" filter. 

We compute the integrals in p.Sp and p.6p using a Gauss-Legendre spectral integrator, and we perform the 
maximization by a Fortran implementation of the Nclder-Mead downhill simplex method [35| . In actual searches, the 
set of all templates forms a grid that should cover the TV-dimensional parameter space. The point on this grid for 
which the FF is maximum singles out the template that best matches the actual waveform impinging on the detector. 
If the FF is above some pre-determined threshold, eg. FF > 0.965, we say that we have a detection. An estimate of 
the number of templates required to detect multi-mode waveforms is given in Appendix [XI 

The detectability of a multi-mode signal by a single-mode template depends, among other things, on the relative 
amplitude of the subdominant modes in the "real" signal. In our simplified signal p.ip we denoted this relative 
amplitude by A. Below we discuss how recent numerical simulations of binary black hole mergers can be used to 
provide such an estimate. 



Estimates of the relative mode amplitude from numerical relativity 

Recent results from numerical simulations of the merger of non-spinning, unequal mass black hole binaries support 
the expectation that the I = m = 2 mode should dominate the ringdown signal However, they also reveal that 
other modes (especially I — m — 3 and I = m = 4) are significantly excited. Whenever the mass ratio q > 2, the 
I = m = 3 mode is excited to roughly one third of the amplitude of the I = m = 2 mode. We expect this estimate of the 
relative excitation to be conservative: perturbative results and numerical simulations indicate that higher multipoles 
should be more excited as the the binary's mass ratio grows, or when the black holes in the binary are spinning. 

Merger simulations show that the projection of the Weyl scalar onto spin-weighted spherical harmonics 
-2Yi„i{0 , (f>) has a circular polarization pattern (see Appendix D of Ref. [l3|)- In the ringdown regime the Weyl 
scalar \I'4 can further be decomposed as a QNM sum, and to a good approximation we can write: 

*4 = J ^ ijlAi\^\e-'/^ _2YirniO , 0) X [siu (x - |m|0) -f I sin (x - \m\cj) + n/2)] . (3.8) 

Ini 
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For simplicity, in this expansion we are considering only the least damped, fundamental mode for each {I , m) pair 
We also approximate spin- weig hted spheroidal harmonics by spin- weighted spherical harmonics, which introduces an 
error of order ~ 1% or less [3a]. The variable x = + '■Po, with ipo a constant, and it is implicitly assumed that all 
frequencies and damping times depend on the angular numbers {I , m). Using the equatorial symmetry of the system, 
the waveform can be rewritten as 



*4 = ; I] c^Milml e~*^' X sin (x ~ mcb) + zl^^^ sin (x - to0 -f n/2)] , (3.9) 

; ,m>0 

where we have defined the following useful quantities: 

= -2Yi,niO,O) + i-iy^2Yi-„,i0,O), (3.10a) 

= -2Ylrn{O,O)-{-iy-2Yl-r,r{0,O). (3.10b) 

Recalling that 'i' 4 — h-^- — ihx and using the large-Q limit, which is usually a good approximation Q, we get 

h+ = -- J2 ^i|m|e-*/^r,+ sin(x-m</>) , (3.11a) 



r 

I ,m>Q 



- E Anrn\e-'/^Y,^, sin ix - mcl> + 7T /2) . (3.11b) 



r 

I ,m>0 



The gravitational wave strain at the detector is given by 

h = h+F+{esAs, V's) + hxFx {OsAs, ^s) , (3.12) 

where ^+,x {Qsi i'Si V's) ^re pattern functions depending on the orientation of the detector and on the direction of the 
source, whose expressions can be found (for example) in Ref. [s^ . In the following, to simplify the notation, we will 
drop the functional dependence oiYj^^ and -f+.x on the angles. 

Let us consider, for simplicity, a two-mode situation. Guided by numerical results for the merger phase (see below) 
we can assume that the dominant modes are I — m — 2 and I = m = 3. Then we get 

h{t) = ^ |e-*/^^^ [sin {ulh + (^22 - 20) Y+F+ + sin (cjf t + (^22 - 20 + 7r/2) r22^^x] - 



A: 



r 

e"*/^^^ [sin {ujfi -I- V333 - 30) Y^j^F+ + sin (ujft + (^33 - 30 7r/2) Y^^Fx\ \ . (3.13) 



^33 
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The relative magnitude of different multipolar components depends on the factor A = Azz / A22 (discussed below) , 
but it is also a function of the angles {9,(1)): that is, it depends on the orientation of the black hole's spin. Let us 
discuss the influence of these two factors, in turn. 

TABLE II: Relative amplitudes of different multipoles, during the ringdown phase. 
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To start with, in Table \U\ we show two different estimates of the relative amplitudes ^33/^122 and ^44/^22 from 
numerical simulations of unequal-mass black hole binaries. The first estimate uses the energy maximized orthogonal 
projection (EMOP) criterion. The idea is to slide a ringdown template of the form (|3.3[) along the numerical waveform, 
and define the starting time of the ringdown as the time maximizing the "energy content" of the ringdown wave, which 
can be defined using a suitable scalar product. Given this starting time, we can also compute the energy radiated in 
a given multipolar component of the ringdown waveform according to the EMOP criterion, i;,™"^ (see Ref. [13 for 
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details). The radiated energy Ef^'-'^ cx (-4™°^)^ QimOJim, where Af^'~'^ is the amphtude of the physical waveform 
h, ujim the vibration frequency of the fundamental mode and Qim its quality factor [5|. Then our EMOP estimate of 
the relative amplitude of different multipolar components is 



/lEMOP 
/lEMOP 



ElmQl'm'^l'T. 
El 'in' Qlm^lr. 



1/2 

(3.14) 



To bracket uncertainties, a second estimate can be obtained by simply computing the ratio of the moduli of 
the waveforms | /if^^'^ | / 1 ^/^^I^ | , where a superscript "peak" means that we evaluate the modulus of the waveform's 
amplitude at the maximum (see Ref. [l3|). We call this the "peak estimate". 

The amplitude ratios depend on the modes being considered and on the binary's mass ratio q, and their functional 
dependence on q can be deduced from the leading-order Post-Newtonian quasicircular approximation discussed in 
Ref. We find that the data in Table HIl are well approximated by the following fitting relations: 

4^ ^ h{l-l/q), (3.15a) 

^22 



-444 , ^ , q^ 

^22 (1+9)' 



~ fc2 + fc3_Ji_^ (3.15b) 



where the values of the fitting constants depend on the estimation method we use. If we use the EMOP criterion we 
get yfcf^OP = 0.303, fcf = -0.0134, fef^^^P = 0.1400. If instead we use the peak estimate, the fitting coefficients 
are fcf"^ = 0.431, fcf"'' = -0.0670, fcf ""^ = 0.2843. 




FIG. 3: ^-dependent angular functions in Eq. (jSTTSj. 

Let us now turn to the angular dependence of different multipolar components. From Eq. (|3.13p we see that the 
relative amplitude of the modes is a complicated function of the sky position and orientation of the source, depending 
on products of the angular functions ^+,x and Y^^^ . A possible way to determine the influence of these angles on 
the relative amplitude of the modes would be to perform Monte Carlo simulations, assuming random distributions 
for the angles. This is beyond the scope of this paper. Some insight can be obtained by plotting the 0-dependent 
combinations appearing in Eq. (|3.13p (Fig. [3]). For simplicity, let us consider the following two cases: 

(i) The angles {ds,4's,4's) are such that — 0. Then the relative amplitude A2/A1 depends on the product 
of two factors: (^433/^22) x (^33/^22)- t^ii^ case, from Fig. [3] we see that the "plus" component of the 
subdominant (/ = m = 3) multipole is never significantly enhanced with respect to the dominant multipole, 
because 1^3^/5^22! ^ 1 ^'^^ values of 9. 

(ii) The angles (9s,4>s,'4^s) are such that F+ = 0. In this case the angular factor (^53/^22) can blow up at the 
equator 9 = 7r/2, so that the "subdominant" I — m — 3 component will actually dominate the waveform for an 
observer located in this direction. 
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In conclusion, the relative amplitude A = A2/A1 is at most '--^ 1/3 for non-spinning mergers with moderate mass 
ratios. Subdominant components can be amplified by the angular dependence of the radiation, but the likelihood of 
such an amplification should be quantified by a more detailed analysis. If the ratio of angular functions is of order 
unity, a relative amplitude A ~ 0.3 can be thought of as a conservative estimate. From both point particle results and 
present-day simulations of spinning binaries in numerical relativity, we expect higher multipoles to be more excited for 
initially spinning black holes and large mass ratios. Below we illustrate the typical effect of subdominant multipoles 
on matched-filtering assuming A — 0.3. 



Fitting factors and event loss due to single-mode templates 

In our calculations of the FF we will assume that the angular dependence of the waveform (|3.13p is such that the 
signal can be simplified to the form of Eq. (j3.ip . This assumption is not general enough. Our discussion above shows 
that it will only be valid for specific source locations in the sky, or for specific orientations of the detector: for example, 
the signal simplifies to Eq. (|3.ip when = or i^x =0, as long as we consider an observer located at some constant 
9 and we appropriately define the azimuthal angles {(pi ,02)- A good strategy to address the general case could make 
use of Monte Carlo simulations. However, the simplified waveform (|3.ip captures many of the important features we 
wish to address in this paper. Based on our discussion of the relative multipolar excitation we assume that mode "1" 
is the fundamental QNM with I = m — 2, and that mode "2" is the fundamental QNM with I = m = 3. We set the 
relative amplitude of the two modes A = 0.3, and to compute the QNM frequencies and quality factors we consider 
a dimensionless Kerr parameter j = 0.6 (a typical value for the end-product of unequal-mass, non-spinning binary 
black hole mergers). 

Fig. [H shows Owen's "minimal match" (1 — FF) [23| for this choice of parameters. The calculation is performed 
for different Earth-based detectors (LIGO, Virgo, Advanced LIGO and EGO) and the minimal match is computed 
as a function of the black hole's mass in the source rest frame, Mq. Thick lines assume that (j>i — 4'2 — 0: roughly 
speaking, this means that the subdominant mode is "in phase" with the dominant multipole. Thin lines assume 
— and (/)2 = TT, a rough way to simulate "dephased" signals. Intuitively we would expect a single-mode template 
to be a worse match for dephased signals. This expectation is confirmed by the fact that the minimal match is usually 
larger when (^2 = (see below for a more detailed analysis). Values above the horizontal dashed line, corresponding 
to a FF — 0.965, correspond to an event loss larger than 10%. This illustrative calculation shows that, for the typical 
relative amplitudes expected from numerical relativity, single-mode templates can produce a significant event loss in 
Earth-based detectors when the black hole mass M > 10^ -^o- This event loss can be fatal for detection of signals 
from very large mass IMBHs (Afo ^ 200 Mq), especially when the I = m = 2 and I = m = Z components are not in 
phase. A more detailed discussion of the dependence of the FF on the phases (0i , ^2) can be found below. 




FIG. 4: Quality factor (top panel) and frequency (middle panel) maximizing the FF, and corresponding event loss (bottom 
panel), for Earth-based detectors. Solid and dashed lines are the template's frequency and quality factor maximizing the 
FF, and the corresponding event loss. Horizontal lines show the frequency and quality factor of the fundamental mode with 
; = m = 2. 
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A low value of the FF is usually accompanied by a large bias in parameter estimation. This is illustrated quan- 
titatively in Figs. [3] and [5l where we plot the quality factor and (dimensionless) frequency maximizing the FF, and 
the corresponding event loss, for different detectors. In Fig. [4] we consider Earth-based detectors of first and second 
generation. To be conservative, we perform the calculation in the "optimistic" case when the I — m = 2 and Z = m = 3 
modes start in phase = (j)2 — 0, A = 0.3). 

The QNM frequency of the dominant mode (horizontal dash-dotted line) is usually estimated with very good 
accuracy by a single-mode filter, except for very large values of the mass (M > 500 Mq for LIGO and Virgo, and 
M > 10"^ Mq for second-generation detectors). Unfortunately, even when the FF is very high the single-mode filter 
produces a large bias in the quality factor of the dominant mode. For all Earth-based detectors we consider, this error 
is > 20% even for low values of the mass {Mq < 100 Mq), when the event loss is quite low and the filter works well 
for the purpose of detection. The bias on the quality factor gets even worse when we allow for a possible dephasing 
of the subdominant multipole (see Fig. [5] below). Being dimensionless, the quality factor of a QNM depends only on 
the dimensionless angular momentum j of the black hole [5,]. Therefore, a large bias in the quality factor means that 
single-mode templates cannot be used for accurate measurements of the black hole's angular momentum, contrary to 
early claims in the literature (see e.g. Ref. [^). Of course, this does not mean that such measurements are not possible. 
Single- mode templates are useful for detection (at least for small black hole masses), but multi-mode templates will 
be necessary if we want to perform precision measurements of a black hole's properties using ringdown waves. 



O 2 





FIG. 5: Same as Fig. |4]for LISA. In the left panel we choose 
also consider a "dephased" QNM superposition with (j>i — 0, ( 



bi = 4>2 = 0, as for Earth-based detectors. 
'2=7!" (green, dot-dashed lines). 



In the right panel we 



Fig. El shows that these remarks remain valid even for LISA, when the SNR is much larger and the detectable QNM 
frequencies much lower. In the plot we consider a source at Dl = 1 Gpc [z ^ 0.2), for which the SNR can be > 10'^ 
0. The left panel shows the "optimistic" case when the first and second QNM signals are in phase. Even in this 
optimal situation the event loss can be as large as ^ 15% for masses Mq ^ 5 x 10^ M©, roughly the measured mass of 
the SMBH at the center of our own Galaxy. When the two QNM signals are dephased (green, dot-dashed line in the 
right panel) the event loss can be larger then 60% for masses Mq ~ 10^ Mq. Notice also that for Mq < 10^ Mq, when 
the single-mode template works well for detection purposes, the bias on frequency and quality factor has opposite sign 
depending on whether the signals are in phase (0i = (/>2 = 0) or dephased (0i = Q, 4'2 = tt). This is no coincidence, 
as we will demonstrate below by studying the dependence of the bias on the phase angles (</>! , (/>2)- 
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Effect of the relative phase of the modes on detection and parameter estimation 

So far we computed the FF assuming (pi — (so that the dominant QNM has maximum amphtude at t — 0). We 
only explored the effect of the relative QNM phase by changing the sign^ of A. In practice the situation is not so 
simple, since the phase of the two (or more) components of the "exact" signal is not known in advance. This problem 
is reminiscent of the analogous problem occurring in matched-filtering detection of inspiral signals (see Appendix B 
of Ref. [111). For both inspiral and ringdown, maximizing the FF (13. 5p over all parameters yields a "best possible 
overlap" which is somewhat optimistic, and therefore not too useful as a detectability criterion. More realistically, we 
should take into account our ignorance of the phase, or phases, of the "true" signal. This can be done by computing a 
"minimax" overlap |3S ]: first maximize the FF (|3.5p over all parameters {A} of the template T{X}, and then minimize 
over the unknown phases of the actual inspiral or ringdown signal h{t). 

To discuss the difference between "best" and "minimax" overlaps in terms of detection and parameter estimation, 
here we perform FF calculations in the (0i , (f)2) plane. We consider, for illustration, two cases: 

(i) An IMBH with Mq = 100 M© as observed by LIGO; 

(ii) An IMBH with Mo = 200 Mq as observed by Advanced LIGO. 

The FF for 0i = 02 = in these two cases is marked by a circle and a square, respectively, in Fig. [T] As in the rest 
of this Section, our signal will be given by the two-mode waveform (|3.ip with A — 0.3 and j = 0.6. This simple model 
is sufficient for our present purpose. A more detailed analysis (taking into account details of the angular dependence 
of the radiation, including a better model of the ringdown signal for spinning mergers, and possibly using waveforms 
from numerical relativity) is a topic for future work. 



Fitting Factor, M=100 M , LIGO 

^ sun 




-3-2-1 1 2 3 

•^1 



FIG. 6: FF as a function of the phase angles of the signal in case (i). 

Fig.[S]shows a contour plot of the FF for case (i). The FF is always larger than 0.965, or equivalently, the event loss 
is always smaller than 10%. For this particular system (and, presumably, for most low-mass black hole ringdowns) 
our ignorance of the phases does not sensibly reduce our chances of detecting the signal. 

The plot shows interesting features, some of which are easy to understand. Under a simultaneous replacement 
(</>! ,4'2) + (2?^ + l)7r,(?!i2 + (2m -I- l)7r), where {n,m) are integers, the FF is unchanged. This is a trivial 

consequence of the fact that the overall sign of the signal p.ip does not affect calculations of the FF. The FF has 
maxima and minima as a function of the two phases. Modulo periodicity, we see that the minimum occurs when 
(j)i ~ 7r/4 and 4>2 ~ —tt/A. This is reasonable: for these values of the phases the signals have comparable magnitude 
and opposite phase at t = 0, so that the I = m = 3 multipole almost exactly cancels out the dominant I = m = 2 



^ For fjii = 0, a sign change in A. is obviously equivalent to setting (p2 = 2n7r (what we referred to as the QNM signals being "in phase") 
or <p2 = (2n -I- 1)tt ("dephased" signals), with n an integer. 
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multipole in the initial (and stronger) part of the signaL This destructive interference produces a signal that is sensibly 
different from a damped sinusoid, reducing the performance of a simple single-mode filter. 



Fitting Factor, M=200 M , Advanced LIGO Matched filtering SNR, M=200 M , Advanced LIGO 

^ sun ^ sun 




FIG. 7: FF (left) and matched-filtering SNR (right) as a function of the phase angles of the signal in case (ii). 

Fig. [7] shows contour plots of the FF and of the matched-filtering SNR pmf, computed according to Eq. (|3.6[) . for 
case (ii). Now the FF is smaller than 0.965 in roughly half of the {(f>i ,(f>2) plane. The event loss ranges from ~ 6% 
to ^ 22%, being larger than 10% in roughly half of the parameter space. From Fig. [T]we can expect that the event 
loss would have been even larger if we had chosen larger values of the black hole mass. The matched-filtering SNR 
was computed assuming that the overall amplitude of the signal Ai corresponds to a ringdown efficiency e^d = 3%, 
and that the luminosity distance Dl = 100 Mpc. Both the FF and the SNR show the (by now familiar) 7r-periodic 
pattern as a function of the phase angles. 



M(a„ M=200 M , Advanced LIGO Quality Factor, M=200 M , Advanced LIGO 

R sun ^ sun 




FIG. 8: Dimensionless frequency (left) and quality factor (right) estimated by a single-mode filter in case (ii). The "true" 
frequencies and quality factors for a Kerr parameter j = 0.6 are: Muri = 0.4940, Qi = 2.9490, MLUii2 = 0.7862, Q2 = 4.5507. 

Suppose that the event loss is moderately large but not so large to prevent a detection, as in case (ii). Then we 
may ask the question: what is the bias in measured parameters when , ^2) maximize the probability of detection? 
In other words: when the template's frequency and quality factor maximize the FF, do they also correspond to the 
"true" frequency and quality factor of the I = m = 2 fundamental mode? Unfortunately, the answer is no. 
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Fig. [8] shows the estimated dimensionless frequency (left) and quahty factor (right) as functions of the phase angles. 
The estimated frequency has relatively small bias, and it always corresponds to the least-damped mode in the pair. 
Results are more interesting for the quality factor. For our chosen value of the Kerr parameter {j = 0.6) the quality 
factor of the I = m = 2 and I = m — 3 modes are Qi = 2.9490 and Q2 — 4.5507, respectively. Comparing with 
Fig. [6l we see that relative minima in pmf (and in the FF) occur when the quality factor "best fits" the I — m = 3 
mode. This is a rather remarkable result: the filter corresponding to the minimax overlap has a quality factor that 
"best fits" the subdominant mode in the pair. Unfortunately, maxima in pmf do not correspond to the filter being 
optimally adapted to the I = m = 2 mode. As the filter tries to maximize the overlap (and the SNR) the estimated 
value of the quality factor becomes significantly biased, and it deviates quite sensibly from the value expected from 
the dominant {I — m — 2) mode. The bottom line is, once again, that single-mode filters can be useful for detection, 
but a multi-mode post-processing will be necessary for accurate measurements of the black hole's angular momentum. 

Another argument in favor of a two-stage search strategy comes from estimating the number of multi-mode templates 
that would be necessary for a detection. Searching for a larger number of modes implies a larger number of free 
parameters, and a correspondingly larger bank of filters. Matched-filtering requires that one covers the possible 
parameter space with a sufficiently fine template mesh, so that the best-matching template lies close enough to the 
true waveform. The distance between ternplates can be quantified in terms of a metric in template space, that was 
introduced by Owen following Refs. [l^, H^. If the mesh is too fine a very large number of templates may be 
required, a computationally expensive option. On the other hand, if the mesh is very coarse the template lying 
"closer" to the true waveform may produce a large event loss. A multi-mode search increases the number of templates 
needed for the mesh to be sufficiently fine. For a single-mode search, Creighton [4^ estimated that ^ 500 templates 
are necessary to reduce the event loss below ~ 10% for LIGO and VIRGO. The same estimate can be used to show 
that ~ 1000 templates would be necessary for single- mode templates with LISA [4l|. In Appendix El we estimate 
that the number of filters Af required for a two-mode template search would be much larger: 



where 5 is a factor of order unity which depends on the detector's frequency span: with our choice of fg we get 
b = 8.3,2.2,1.6,1.2^^1.6 for LISA, EGO, Advanced LIGO, LIGO and Virgo, respectively. Using better template 
placing techniques [25| or imposing constraints on the functional form of the QNM frequencies (Appendix |^ could 
help reduce computational requirements. A two-stage search seems to be a good compromise between performance 
and computational costs. Single-mode templates can be used for detection. Given a detection, multi-mode templates 
or Prony methods [l^l should be used for parameter estimation. A larger number of templates also means that the 
threshold for detection must be set higher, because there is a larger false alarm probability. Hierarchical searches 
or other techniques could play an important role in this regard |42l |. In any case, such a large number of templates 
may not be a problem by the time advanced detectors are in operation. Improved computer performance and the 
use of large-scale computational projects, such as Einstein@IIome [435, could be sufficient to overcome computational 
difficulties within the next decade. 



So far we looked at the event loss and bias in parameter estimation due to the use of single-mode templates to 
detect multi-mode signals. The discussion assumed that the gravitational wave signal is composed of at least two 
QNMs having roughly comparable amplitude. The question we address here is the following: how can we tell if 
there really are two or more modes in the signal, and can we resolve their parameters? If the noise is large and the 
amplitude of the weaker signal is very low, or the two signals have almost identical frequencies^^ the two modes could 
be difficult to resolve. This issue is particularly significant since no-hair tests using ringdown 5] require the presence 
and resolvability of two or more modes. Roughly speaking, the first mode is used to measure M and j by inverting 
the relations /i = fi{AI ,j) ,Qi = Qi{M ,j) (see Appendix El for more details); the second mode can then be used 
to test consistency with the Kerr geometry. 

Here we really address two different issues. We first assume that there are indeed two modes in the signal, and 
we discuss criteria to resolve their frequencies and damping times. This discussion parallels that in Ref . ,51] , updating 
estimates of the relative mode excitation on the basis of recent results from numerical relativity, and correcting a typo 
in that paper. Then we introduce a rigorous criterion to resolve amplitudes: that is, we compute the minimum SNR 
such that one can decide by the presence of two modes in a given ringdown signal. 

Resolving frequencies and damping times. A crude lower limit on the SNR required to resolve frequencies and 
damping times was presented in Ref. Q. The analysis uses the statistical uncertainty in the determination of each 




(3.16) 



IV. 
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frequency and damping time, which a standard Fisher Matrix calculation^ estimates to be @ 

1/2 



^ / /f(3 + 16Qt) 
2 / (3 + 4Q2] 



AihQi 



/i(l + 40?) /2(1 + 4Q2)_ 

1/2 



/i(l + 4Q?) /2(1 + 4Q2) 



(4.1a) 



(4.1b) 



These errors refer to mode "1" in a pair. By considering the "symmetric" case (/"i = 02 = 0, the errors on /2 and T2 
are simply obtained by exchanging indices (1 ^ 2). The expression above holds in both the FH and EF conventions, 
assuming white-noise for the detector, but modes "1" and "2" must correspond to different values of I or m (in the 
nomenclature used in Ref. 0, the QNMs must be quasi-orthonormal) . 

A natural criterion (d la Rayleigh) to resolve frequencies and damping times is 

I /i - /2 1 > max(CT/j , CT/J , In - T2 1 > max(cr^j ,(7^^). (4.2) 

In interferometry this would mean that two objects are (barely) resolvable if "the maximum of the diffraction pattern 
of object 1 is located at the minimum of the diffraction pattern of object 2" . We can introduce two "critical" SNRs 
required to resolve frequencies and damping times. 



/ _ max(po-/, , pcT/J ^ _ max(pg^,,pg^J 

^'crit ~ \ r r \ I Pcrit ^ i _ i ' \^-'^) 

1/1 -/2I \Ti-T2\ 

and recast our resolvability conditions as 



p > pcrit = min(p£.ij,p^,.it) , (4.4a) 
p > pboth = max(p4it,p^,,it) . (4.4b) 

The first condition implies resolvability of either the frequency or the damping time, the second implies resolvability 
of both. 

Resolving amplitudes. A related question is: how large a SNR do we need to confidently say that we have detected 
a multi-mode signal, and to resolve two signals of different amplitudes? Suppose again, for simplicity, that the true 
signal is a two-mode superposition. Then we expect the weaker signal to be hard to resolve if its amplitude is low 
and the detector's noise is large. 

In Appendix |B] we quantify this statement by deriving a critical SNR for amplitude resolvability based on the 
generalized likelihood ratio test (GLRT), Pglrt- The derivation of this critical SNR, which is given in Eq. (|B12|) . 
is based on the following simplifying assumptions: (i) using other criteria, we have already decided for the presence 
of one ringdown signal, and (ii) the parameters of the ringdown signal (frequencies and damping times), as well as 
the amplitude of the dominant mode, are known. In practice the latter assumption is not valid. For this reason, our 
estimates of the minimum SNR needed to detect more than one mode should be considered optimistic. 

Fig. [5] compares the critical SNR pglkt, as defined in Eq. (|B12p . and the two different criteria for frequency 
resolvability of Eq. (|4.4p . All quantities are computed as functions of the binary's mass ratio q. The angular momentum 
j of the final black hole is computed using the fitting formula derived in Ref. for the angular momentum of the 
black hole resulting from unequal-mass, nonspinning binary black hole mergers: j = 3.352?7 - 2.46177^, where the 
symmetric mass ratio t] = q/{l + q)^ . This value of j is then used to read QNM frequencies from numerical tables. 
We assume the dominant mode to be the fundamental I = m = 2 QNM and for the subdominant mode we take the 
fundamental QNMs with l = m = 3 or l^m = 4. To compute the relative amplitude A{q) for different mass ratios 
we use the EMOP estimate of Eq. ()3.15a|) . 

The plot shows that pcrit < Pglrt < Pboth for all values of q. Therefore, given a detection, the most important 
criterion to determine whether we can carry out no-hair tests seems to be the GLRT criterion. If p > pglrt we can 
decide for the presence of a second mode in the signal. Whenever the second mode is present, we also have p > Pcnt- 
that is, we can resolve at least the frequencies (if not also the damping times) of the two modes. A SNR p 30 — 40 
is typically enough to perform the GLRT test on the Z = to = 3 mode, as long as q > 1.5 or so (equal-mass mergers 
should be quite rare anyway). By looking at Fig. [2] we conclude that not only LISA, but also advanced Earth-based 
detectors (Advanced LIGO and EGO) have the potential to identify Kerr black holes as the vacuum solutions of 
Einstein's general relativity. 



^ Our Eq. I l4.1al l corrects a missing factor of 2it in Ref. 
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FIG. 9: Minimum SNR required to resolve two modes, as function of the binary's mass ratio q. If p > pglrt we can tell 
the presence of a second mode in the waveform, if p > pcrit we can resolve either the frequency or the damping time, and if 
P > Pboth we can resolve both. Mode "1" is assumed to be the fundamental mode with / = m = 2; mode "2" is either the 
fundamental mode with I = m = ?> (solid lines) or the fundamental mode with Z = m = 4 (dashed lines) . 

V. CONCLUSIONS 

In this paper we analyze the detectability of ringdown waves by Earth-based interferometers. Confirming and 
extending previous analyses, we show that Advanced LIGO and EGO could detect intermediate-mass black holes of 
mass up to ~ 10^ Mq out to a luminosity distance of a few Gpc. 

Using recent results for the multipolar energy distribution from numerical relativity simulations of non-spinning 
binary black hole mergers to estimate the relative amplitude of the dominant multipolar components, we point 
out that the single-mode templates presently used for ringdown searches in the LIGO data stream could produce 
a significant event loss (> 10% in a large interval of black hole masses). A similar event loss should affect also 
next-generation Earth-based detectors, as well as the planned space-based interferometer LISA. 

Single-mode templates are useful for detection of low-mass systems, but they produce large errors in the estimated 
values of the parameters (and especially of the quality factor). We estimate that, unfortunately, more than ~ 10^ 
templates would be needed for a single-stage multi-mode search. For this reason we recommend a "two stage" search 
to save on computational costs: a single-mode template could be used to detect the signal, and a multi-mode template 
(or even better, Prony methods [13]) could be used to estimate parameters once a detection has been made. 

In Appendix [B] wc introduce a criterion to decide for the presence of more than one mode in a ringdown signal. 
By updating estimates of the critical signal-to-noise ratio required to resolve the frequencies of different QNMs using 
results from numerical relativity, we show that second-generation Earth-based detectors and LISA both have the 
potential to perform tests of the Kerr nature of astrophysical black holes. 

In the future we plan to use numerical waveforms (possibly including spin effects) to refine our estimates. We also 
plan to carry out Monte Carlo simulations to study the information that can be extracted on the source position and 
orientation using a network of Earth-based detectors. The possibility to constrain the black hole spin's direction from 
the multipolar distribution of the merger-ringdown radiation should be particularly interesting (eg. for coincident 
electromagnetic observations of jets that could be emitted along the black hole spin's axis). 
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APPENDIX A: NUMBER OF TEMPLATES FOR A TWO-MODE SEARCH 



Two-mode ringdown waveforms of the form (|3.ip depend of five parameters: two quality factors Qi,Q2, two 
frequencies /i , /2 and one relative amplitude A between the different modes (for simplicity, here we set 0i = '/>2 = 0) . 
For matched-filtering with ringdown waveforms of unknown frequency and quality factor we must lay down a "mesh" 
covering the parameter space with some pre-defined precision (that can be translated to a pre-defined minimum loss 
of SNR) [U, [3^, ll^] ■ We follow Owen's formalism to estimate the necessary number of templates. The distance 
between nearby templates, which defines the mismatch between the filters, can be computed in terms of the metric 

ds^ ^ g^^dx'^dx^ = (gi,Q2,/i,/2,-4), (Al) 

where (for large or moderate values of Qi and Q2) the metric coefficients are well approximated by 



^ /2 (/2Q1 + 2A'hQ2) ^ A^h {2f2Qi+A'hQ2) 

^'^"'^^ 8Qi(/2Ql+^VlQ2)" 8Q2if2Ql+A^hQ2f ' 

Sf^f- fim{hQi+Ay,Q2)' f^'^Hf2Qi + A^hQ2)' ^ ' 

_ flf2QlQ2 _ A^hf2 . 

" 2 ihQ, + A^f.Q^f ' - 8 if2Q^ + ^ViQ.)^ ' ^ 

^ /2 (/2Q1 + 2A'hQ2) ^ A^hQ2 .^2dl 

^"^'^^ 8/i(/2Qi+^ViQ2)" ^"^'^^ 8{f2Qi + A^hQ2f' 

^ A^f2Qi A^h{2hQi+A^hQ2) 

-4/1/2Q2 _ ^Ahf2Qi_ 

4 (/2(3i +-4^/1(32) 4:[f2Qi+A^fiQ2) 

A^QiQ2 Af2QiQ2 . 

''''' HhQ.+Ay,Q2f' '''^ Hf2Qi+A^hQ2f' ^ 



AfiQiQ2 
4{f2Qi+A^fiQ2) 



9f2A = ~ w . ^ , ^9 r ^ .2 - (A2h) 



Requiring a loss of no more than 10% in the event rate due to a mismatched template (i.e, the minimal match MM 
[H, [40| should be at least 0.97) we get an estimate for the number of templates we need: 

/ dQ 1 dQ2dfi df2 yd et I |g^, 1 1 ^^^^^^^f 0.03 
32 [(1 - MM)/5]' 



Here 6 is a factor of order unity which depends on the detector's frequency span. We get h = 8.3 , 2.2 , 1.6 , 1.2 , 1.6 for 
LISA, EGO, Advanced LIGO, LIGO and Virgo, respectively. In deriving this number we assume the frequencies to 
be searched for are those of interest for each of the detectors (3 x 10^^ ^ / 1 for LISA, 10 < / < 2000 for EGO, 
20 < / < 2000 for Advanced LIGO and Virgo, 40 < / < 2000 for LIGO), that the quahty factor varies between 
and 20 for all modes likely to be detected f5|, and that the relative amplitude A varies between 0.01 and 100. Our 
estimates are not strongly sensitive to the relative amplitude: assuming < .4 < 1 yields a total number of templates 
which is roughly half the above number. 

For the single-mode case, setting Q2 — f2 = A = we get the following metric: 

^ 8^^^' ~ i^^^'^^ + ■ ^^^^ 

The number of templates Af ^ 6 Qmaxlog (/max//min) 0- In particular we get A/" ~ 460 for LIGO and Virgo [40| . 
and A/" ~ 1000 for LISA [4l| (a huge improvement in terms of computational requirements). 

The formalism and numbers presented in this section are valid for a general ringdown signal: no constraints were 
imposed on the QNM spectrum. A possible approach to reduce the number of templates is to assume that the 
source is a general relativistic black hole. In this case, we are left with only three intrinsic parameters: the mass and 
angular momentum of the black hole and the relative amplitude between the modes. Alternatively, we can choose the 
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independent parameters to be one quality factor Qi, one frequency /i and the relative amplitude of the modes; the 
quality factor and frequency of the second mode, Q2 and /2, can be thought of as functions of Qi and /i. 

More explicitly, we find that for rotations ^ j ^ 0.98 the frequencies and quality factors of the fundamental mode 
with I = m ^ 2, 3, 4 are well approximated (to within ~ 6% or better) by 

3.2+4.91 

2TTMfm ~ 0.74 + 0.39; - (0.78 + 0.18/) (1-j)^™- ^ (A5) 
Qm ~ 0.26 + 0.22/ +(-0.36 + 0.88/) (l-j)~° '*V (A6) 

By inverting these relations for (say) the / = m = 2 mode, one can infer [j, M) and compute the frequency and 
quality factor of the modes with I — m — 4. More generally, in Ref. 5] the frequencies and quality factors of the 
first three overtones (for / = 2, 3, 4 and all values of m) were fitted by functions of the form 

2^A//i„.„ = /i+/2(l-j)^% (A7) 
Qhnn = qi+q2il-jy'. (A8) 

Here the constants fi and qi depend on {l,m,n) (see Tables VIII-X in ^]), and the fits are accurate to better than 
4%. By using these fits or (more precisely) numerical QNM data we can express any subdominant mode in terms of 
the dominant mode, and possibly reduce the number of templates. 

Another possibility to reduce the number of templates is to restrict the parameter space by using information 
derived from measurements of the inspiral waveforms. For example, if we had a reasonably accurate measurement of 
the masses and spins of the binary members we could significantly restrict the possible values of the mass and spin of 
the final black hole to be searched for. 



APPENDIX B: AMPLITUDE RESOLVABILITY 



The purpose of this Appendix is to estimate the minimum SNR required to test the hypothesis that a second 
mode is present in a ringdown waveform. The derivation is based on the following simplifying assumptions: (i) using 
other criteria, we have already decided for the presence of one (dominant) damped sinusoid in the signal, and (ii) the 
parameters of the ringdown signal (frequencies and damping times), as well as the amplitude of the dominant mode, 
are known. In practice the latter assumption is not valid, so our estimates of the minimum SNR should be considered 
optimistic. 

The question of whether one or two damped sinusoids are present in the signal can be stated in statistical terms, as 
follows''. Let w{t) be a zero-mean gaussian white noise process, and define y{t) = s{t) ~ Aihi{t) to be the difference 



between the actual signal s{t) and the dominant normalized QNM signal hi{t) = e 



-{^fi/Qi)t 



sin(27r/it) (for simplicity. 



in the present discussion we set (jji — (j)2 — 0) . Denote by Tii the hypothesis that the signal contains only one damped 
exponential in noise, and by Ti.2 the hypothesis that the signal consists of two damped exponentials in noise: 



Hi : y{t) = w{t) 

H2 : y{t) =A2 h2{t)+w(t) 



(Bl) 



We write the normalized second QNM as h2{t) = e"(''-''2/'32)* sin(27r/2i). 

We assume that we do not possess any prior information on the possible values of A2, besides the fact that A2 > 0. 
The general structure of composite hypothesis testing is involved when unknown parameters (in this case, A2) appear 
in the probability density function. We follow Milanfar and Shahram [3, HE] and consider the generalized likelihood 
ratio test, which proceeds by computing first the maximum likelihood (ML) estimates of the unknown parameters. 
These estimates are then used to form Neyman-Pearson detectors. We note that this approach has been compared 
very favorably against other standard tests in Ref. [44| . 

Suppose that the signal s(t) is sampled at discrete times tk (k = 1,...,7V) and that the corresponding sample 
values are — s{tk), yk = Sk — Aihi{tk). If we assume additive white noise with variance a, the probability of the 
observed data under the hypothesis of a second damped sinusoid of amplitude A2 in the signal is 



N 



N 



fe=i 



2tt 



(Vk - A2 h2k) 
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2a2 



(B2) 



See also the work by Milanfar and Shahram l44l |45|| . 
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To decide on hypothesis Tii or 7^2, we start by computing the ML estimate of the unknown parameter A2- 



N 



maxpA2 = maxlnp^, = miny^(j/fc - A2 h2kf = min ||Y - A2 H2IP , (B3) 

k—1 

where Y and H2 are column- vectors of the samples yk and /12 k , respectively, a superscript T stands for "transpose" , 
and ||V|p = V'^'V denotes the norm of a vector V. The (unconstrained) ML estimate A2 of the parameter A2 is then 

We wish to test the hypothesis that A2 > against the hypothesis that A2 — 0. There is no general-purpose test 
to do this, but a powerful test is to compute a likelihood ratio and to find the maximum in both the numerator and 
denominator: this is called the generalized likelihood ratio test (GLRT). Some algebra shows that 

j.(Y) ^ In ^^-n2PA2{y.) ^ 1 / H|V \ ' 
ma-x-Hi PA2=o{yx) 2\a\\B.2\\ 



For any given data set Y, we decide on 7^2 if ■\/2r(Y) exceeds a specified threshold 7: 

H?^Y , , 

While it may seem troublesome to use the unconstrained ML estimate to form the GLRT, in fact, due to the assumed 
positivity of A2 , the detector structure is effectively producing a one-sided test and therefore this is in fact a uniformly 
most powerful detector (46j . 

The choice of 7 is motivated by the level of tolerable false-positive rate [H, . The detection rate Pd and false 
alarm rate P/ for this detector are related as 

Pd - QiA2 77 + 7) = Q{A2 V + Q'\P})) , (B7) 

where 

IIH2II 

V = • (B8) 

a 

Here Q denotes the right-tail probability function for a standard Gaussian random variable (zero mean and unit 
variance): 



Qix) 

From (IB7I) we find 



.00 ^ 


r 91 

W 


/ exp 




Jx 


2 



dw . (B9) 



Q'\Pd)-Q'\Pf)^A2V- (BIO) 
This last expression can now be put in a more convenient form, by defining the output SNR as 

^J|^A+^2H2||_ ^^^^^ 
(T 

Using (|B8|) and (jBlOp the relation between the minimum resolvable A — A2/A1 and the required SNR can be made 
explicit. The critical SNR for hypothesis testing pglrt is 

PGLRT = [Q-\Pd) - Q-\Pf)] ""^il+^^l^^" . (B12) 

For the calculations in Fig.lHlwe set Pd — 10^^, Pf — 0.99. Had we chosen Pd = 0.1 and Pf ~ 0.9, the critical p would 
have decreased by a factor ~ 2.6/4.6. For more stringent false alarm rates (say, with Pd = 10^^ and Pf ~ 0.99) the 
critical SNR would have increased by a factor ~ 7.1/4.6. 
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